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SUMMARY 



This investigation was carried out in an effort to 
express in analytical terms the relationship between a 
pilot's Input and output while flying a Eyan Navion 
Aircraft during the performance of a specified maneuver 
in the longitudinal plane. The maneuver chosen for the 
pilot to perform was a simple air speed change in which 
the pilot was required to reduce the aircraft velocity 
from a given value to a lower assigned value using only 
the elevator. He was further required to do this in a 
smooth and expeditious manner and to arrive at his new 
air speed without overshoot or undershoot. The mathematical 
relationship may be expressed as the ratio of his output, 
(the manner in which he positions the elevator), in response 
to his Input (the desire to change air speed a given 

number of miles per hour). This ratio, H » § , is called 

6 

the human transfer function. 

In order to determine the human transfer function, 
it was necessary to divide the investigation into three 
parts. 

In Part I, an expression for the ratio of the aircraft 
velocity response to an elevator forcing function in the 

form of a step was obtained. This ratio, A s u s — , can 

% fc 

be considered the aircraft transfer function, and can be 
expressed with good accuracy by a second order differential 
equation whose characteristic equation yields roots from 



which the period and time to damp to half amplitude of 
the aircraft phugoid mode may be calculated. 

In Part II, several mathematical expressions for the 
human transfer function were obtained. Ir. arriving at all 

of the expressions for, H * & , certain simplifying 

e 

assumptions were made. It was considered that the pilot 
was merely a "black box" In the system, and that his inputs 
were limited to a visual realization of the error (air 
speed change desired) and simple functions of the error. 

No attempt was made to consider any time delay in the 
pilot's resoonse, nor were any neurc-muscular effects 

taken Into account, A differential equation containing 

• • • • « * 

terms of S , 8 , fc , £ , 6. and & , and an integro-dlf ferentlal 
equation containing terms of & , 8 , £ , 8. and dt were 
each solved simultaneously with the aircraft transfer 

functions on an analog computer and time histories of velocity 

* 

and elevator response were recorded. Comparison of these 
responses with those taken In flight showed that the 
differential equation gave slightly better correlation with 
the velocity response than did the integro-dlf ferentlal 
equation. However, this difference was slight. Comparison 
of the computer solution of the elevator response with 
that recorded in the airplane showed that the lntegro- 
differentlal equation gave the better correlation. It Is 
felt that the integro-dlf ferentlal equation represents the 
pilot's transfer function reasonably and well. Recommend- 
ations were made to check the validity of both expressions 



under varying conditions. 

In Part III of the investigation, a mock-up 
simulating the elevator controls of the aircraft was 
constructed. Pore and aft movement of the control column 
caused a positive or negative voltage output from the 
mock-up which corresponded to the elevator angle displace- 
ment. This voltage was used as a forcing function on the 
simplified equations of aircraft motion which were set up 
on the analog computer. The velocity output on the computer 
was displayed on a meter located in front of the pilot. An 
attempt was made to correlate the pilot's movement of the 
elevator in accomplishing speed changes on the computer 
with that occurring in the air and with the results obtained 
In Part II. No satisfactory comparison could be made 

because of the rough design of the mock-up. It is felt 

% 

that with an Improved raock-up design and with certain 
recommended changes In the test procedure a fruitful com- 
parison might be made. 



3Y11B0LS AND SIGH CONVEIITIOH 



£ e Absolute elevator angle in degrees. Positive down. 

6 Perturbation elevator angle in degrees taken from 
an arbitrary reference, Dov;n elevator positive. 

& d^/dt and d <4 /dt^ respectively. 

£ The error in mph from a desired airspeed reference. 
Positive error if the speed is higher than desired. 

<£, i d£/dt and d*£ /dt 2 respectively. 

P Aircraft period in seconds. 

7^/2 Time to damp to half amplitude in seconds* 

D Tire operator d( )/dt. 

u Hon dimensional velocity change change, AY/ V# 

V Velocity in fps or mph as indicated. Subcripts "i” 

and ”f n are initial and final velocities respectively. 




Stick force in pounds, 

= Partial derivative of the moment coef icient 

with respect to angle of attack. 



Slope of the lift curve. 

c mdO “Jtz - P art i a l derivative of the moment coefficient 
^ with respect to 7“-e-. 




W -i*c 



Partial derivative of the drag coefficient 
with respect to angle of attack. 

Partial derivative of the moment coefficient 
with respect to T<c . 



C m ^ Elevator power. 
oC Angle of attack. 

d€/d«c hate of ch.ange of doxmwash angle with angle 
of attack. 



m 



I lass of the airplane in slugs 



/“• = n/p3c 
T = m/josv 

Iy I'loment of inertia about the Y-axis in slug ft* 
h = 

A => (f) a i^ cr aft = iiircra ^ Transfer Function 
H = (f-) liuri ~ n = - unan Transfer Function 



NAVION SPECIFICATIONS AND NOTATIONS 
Wing area, S = 184.2 ft^ 

Horizontal tall area, S t = 43 ft^ 

Elevator area - 13.04 ft^ 

Wing aspect ratio, A, = 6.04 
Wing span, b, = 33.38 ft. 

Tall span, b t = 13.17 ft. 

Wing MAC, c, - 5.7 ft. 

Tall MAC, c t , - 3.34 ft. 

Tall length, l t = 15.04 ft. 

Incidence, wing, l w , = 2 © ( a t root), 10 (tip) 

Incidence, Tall, 1 ^., - -30 

Airfoil, Wing, Root: NACA 4415 R, Tip: NACA 6410 R 
Airfoil, "Tail, NACA 0010 
Dihedral r- 7.5° 

Wing Taper ratio, A, s o 54 

Wing aerodynamic center, a.c., .242 

Elevator: 

Root chord =1.5 ft. 

Tip chord = 1.0 ft. 

Stabilizer: 

Root chord =2.5 ft. 

Tip chord "1.67 ft. 

Thrust axis is, parallel to fuselage reference line. 



Cl = lift coefficient. Subscripts a, w, t denote 
airplane, wing and tail, respectively 
C ra - moment coefficient 

C X 

v * tail volume ■ t t * # 6l6 

S w c 

c.g. = center of gravity 

= angle of attack. Subcripts w, t, and p denote 
wing, tail and propeller thrust axis respectively. 
X C g = c.g. position in percent of MAC 
Xac = a.c. position in percent of MAC 
Sp = propeller disk area = 38.5 ft^ 

N 0 - control-fixed neutral point. 



INTRODUCTION 



During the past ten to fifteen years the advancements 
In the engineering and other physical sciences have far 
surpassed those of the biological sciences. As a result, 
new equipment of every sort is being developed, built, 
and put Into use with no scientific assurance that, from 
the standpoint of the human operator, it is completely 
useful or habitable. As a consequence, the field of 
human engineering with its task of increasing the effective- 
ness of a man-machine system by treating it as a unified 
system has become of Increasing Importance. Organizations 
in various universities, colleges, civilian enterprise, 
and the military are attempting to apply the known principle 
of psychology, physiology, anthropology, and medicine, 

at 

as well as engineering techniques, to the man-machine 

-t 

relationship and equipment design. Further, they have 
set up elaborate programs for research in order to increase 
their knowledge in these fields. 

The problem of aircraft control illustrates the 
complexity of the human response. In controlling the air- 
craft, the pilot receives certain data from his flight 
instruments which he judges, compares, accepts or rejects, 
and then moves his controls accordingly. This control 
movement changes the information sent to the pilot by his 
instruments, and the pilot must continually re-evaluate 
this information in an effort to make the aircraft perform 



in the manner which he desires. In actuality, this is 
more than a simple, error-sensing, closed loop servo 
system, for the pilot is capable of fulfilling (and does) 
the additional functions of examining the information 
received from his instruments in the light of their rates 
of change and their time duration. This means that he 
is acting in some respects as a computer which senses errors 
and their derivatives or integrals. Theoretically it 

i 

should be possible to express the pilots relationship 
between the inputs and outputs of the system (i.e. the error 
in some Instrument from a desired value and the pilot's 
control movement to correct the instrument reading) as 
some mathematical relationship corresponding to a closed 
loop servo system which Includes one or more types of 
feedback. 

A mathematical expression of this sort is tantamount 
to considering the pilot as a "black box" in the system, and 
does not consider the biological, psychological, or physio- 
logical principles which enable him to perform these 
functions. However, such a mathematical expression should 
be of interest in attempting to explain why the human acts 
as he does. 

The authors of this report have elected to use the 
"black box" approach in their investigation of the relation- 
ship of the human and the aircraft. One simple task was 
assigned to the pilot. He was required to decrease his 

air speed from a given value to another specified velocity, 



smoothly and without overshoot. 

In an effort to express the pilot's actions in 
performing this task by an analytical expression, or 
transfer function, the problem was divided into three parts. 

I Determination of the aircraft transfer function. 

II Determination of the human transfer function. 

III Check of the validity of the human transfer function 
by a Mock-Up. 



PART I 



DETERMINATION OP THE AIRCRAFT TRANSFER FUNCTION 



THEORY 



In Refs. 1 and 2 It Is shown that the aircraft 
equations of motion can be combined to give a character- 
istic equation of the form: 

CJ\ + Ci\ + CA + 0 o = O (1) 

Solution of this equation (confirmed by flight 
experience) has shown that the aircraft has two modes of 
longitudinal motion. The large pair of roots correspond 
to a short period, heavily damped mode, while the small 
pair of roots correspond to a lightly damped, long period 
oscillation called the phugoid. The response of an air- 
craft of normal conf iguration to an elevator forcing 
function is such that the short period mode is so heavily 
damped that either the pilot is not aware of its presence, 
or, if aware. Is unable to make a correction before the 
motion is completely damped out. This is particularly 
true of the aircraft velocity response which, in general, 
experiences no visible change. However, the presence of 
the phugoid with its long period (30 to 60 seconds) and 
light damping (time to half amplitude of 30 seconds or 
longer) is easily detected by the pilot as a change of 
airspeed and as a change of attitude. 

Thus, as far as the pilot is concerned, the Important 
part of the aircraft velocity response to an elevator 



forcing function is the phugold mode. The characteristic 
equation, when the short period mode is neglected, then 



reduces to the form: 

K £ ' 

'Z 



K X +• V<,„ — O 



( 2 ) 



If we neglect the short period mode and assume that 
the aircraft has a transfer function of the form: 

U - ' 

S 



( cx T)^ + bD t vO) 



(3) 



£> - elevator 

angle applied, degrees 



- AV 
where u. ~ \j 

_ _ \J - initial velocity, 

^ mph 

Then the characteristic equation would be: 



CS. 



* 






Yx = o 



(4) 



which is of the same form as Equation (2) 
The roots of Equation (*0 would be: 






2.cx 



- 



+ 



o. 



( 5 ) 



If the aircraft were very lightly damped, then, 

-te- «. ^ 

2o, a. 

and the equation would reduce to: 



K* = 



2cx - M 



v\ 

CX 



( 6 ) 



The expressions appearing in Equation (6) are related 
to the period and time to damp to one half amplitude 
by the following equations (Ref. 1): 




VJO 









Sec. . 



(7) 

(8) 



If a known elevator forcing function were applied 
to the aircraft in flight, and the time history of the 
aircraft velocity could be recorded, then the period and 
the time to damp to one half amplitude could be measured 
from the time history of the velocity. From this data, 
and knowing that 





steady state or final 



a 



steady state or final 



the value of the unknown coefficients of the aircraft 
transfer function, Equation (3), could be computed from 
Equations (7) and (8). 



PHOCEDUBE AND DISCUSSION 



As indicated in Appendix A, the aircraft was instru- 
mented to record velocity and elevator angle. In order 
for the pilot (seated at the left hand controls) to apply 
an elevator step function to the aircraft in flight, a 
positive stop to control movement was erected on the 
right hand controls, A sketch of this stop is shown below. 



The stop was adjusted so that the stick movement corre- 
sponded roughly to an elevator change of -.7 degrees (up 
elevator). This elevator deflection produced an airspeed 
change of approximately 20 miles per hour. 

The test flights were conducted in the following 
manner. All instrumentation was calibrated prior to take 
off, and upon return. Plights were made during the very 
early morning in order to have smooth air condldtlons. 
Each run was commenced from a base pressure altitude of 







Posa WF- 
S^OP 






4800 feet in order to give an average pressure altitude 
of 5000 feet during the run. Temperature was recorded 
at 5000 feet in order to correct the recorded calibrated 
air speeds to true air speeds. The pilot established a 
base calibrated air speed of approximately 120 mph at 
4800 feet and trimmed the aircraft. After the recording 
mechanism was started and checked, the pilot moved the 
yoke sharply rearward until it met the positive stop, 
and then held it firmly against the stop until all air- 
craft oscillation had ceased and the aircraft was again 
in steady flight at a lower airspeed and a less positive 
elevator angle. 

A series of eight runs of this nature were made and 

the time histories of the aircraft velocity response were 

recorded. As no effort was made to start each run from 
% 

exactly the same air speed, and as the step function could 
not be reproduced exactly from rim to run, the aircraft 
did not settle down to the same speed on each run. 

Examination of the time histories of the elevator 
response showed that there was some random variation of 
the period and time to damp to half amplitude from run 
to run. Hov/ever, all Periods were 36 seconds ^ 1 second, 
and all Times to half amplitude were 2 5 seconds "t 1 second. 
The Via 

s/vV\er e ~ Vxrxoa 

- \ Vn\tyoA ~~ yl f \ rvcxA 
M X r\\ V vv\ 

was computed for each run and all results lay in the range 



4.70t .10. Values of P, T|, and b/u for the eight runs 
are tabulated In Table I. Two typical time histories of 
the velocity response taken from the recording oscillo- 
graph record are plotted as Figs. 2 and 3 bo show the 
phugoid mode. 

After examination of the results in Table I, it 
was elected to choose the values obtained in Bun 2172 for 
computing the aircraft transfer function. These values 
are reproduced below for convenience: 



= 132 mph 

V f = IO 8.5 mph 
AV = 22.5 mph 

8 = . 80 degrees 



u = _V = 22.5/132 = .1705 
Vi 

8 = K « 4.70 
u 

P a 36 seconds 
T| - 25 seconds 



The constants of Equation ( 3 ) were then computed from 
Equations f?) and (8). 



p= Wf? 






ex *=■ 



bV~ ( A >10^1 
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' 0/2.CV. 



T 



N / x. 



IS 



^.50 



The equation (3) could be written: 



m _ 

% 



\ 



V ^>.50 D -V- 



(9) 



(\5Vfe^' 






Expressing the transfer function in terras of £^V, we 
have: 

£>V _ _ L_ — 

* .ObA-A-D -v .0355) 

It is seen that the above equation satisfies the initial 
and final physical conditions: 




at t * O 

kV = O ' *6 = 0 



Owt 



t - ‘b'tcx'te 

£ V - 2."2-.b 



Now if we consider that when flying at the higher 

*c 

airspeed, 130 mph, and at an elevator angle of .80 degrees, 
we desire to- decrease the airspeed to IO 8.5 mph by the 
application of a step function elevator input, then the 
elevator angle must be decreased to zero by the step 
function. Thus if we define 

error - £ - V t - V ^ 

elevator input - S " ^ ~ 




the aircraft transfer function may be written 




where £ represents the difference between 
any speed and the speed desired 
S represents the elevator step function 
required to attain the desired speed. 
This is the form In which the aircraft transfer function 
is used during the remainder of this report. 

This expression was then checked by means of an 
analog computer. The velocity time history was recorded 
on a Brush Recorder In response to an elevator step 
function voltage corresponding to .80 degrees. The period 
and time to half amplitude were found to be seconds 
and 25 seconds, respectively. 



RESULTS 



Flight tests indicated, that the phugold mode was 
relatively constant with a period of 36 seconds and a 
time to half amplitude of 25 seconds. Further, they 
show that the ratio of the elevator step input to the 
speed change resulting from that Input Is constant. 

The aircraft transfer function as determined by 
flight tests and checked by an analog computer may be 
expressed as: 

^= 6 ^= . ( 11 ) 

* ( "D -r ,0fe*4T> -V .O-ISV) 



CONCLUSIONS 



The pilot in flying the aircraft does not recognize 
the presence of the short period mode in the velocity 
response of the aircraft to an elevator input. The 
recorded time history of the velocity response does not 
show the short period mode. 

The pilot is aware of the phugold mode caused by an 
elevator disturbance. 

The transfer function of the velocity of the aircraft 
in response to an elevator forcing function in the form 
of a step can be closely approximated by a second order 
differential equation. 



PART II 

DETERMINATION OP THE 



PILOT'S TRANSFER FUNCTION 



THEORY 



In an effort to express the pilot's actions to 
a desired change in air speed in terms of an output and 
an Input, he may be considered part of a servomechanism. 
This system is a simple closed loop system, sketched belov/: 



tr 




P,/ot 


& 


A i rcrdft 


*5 


? 


H 




A 



V ln is a step function Input of velocity (or the 
desired velocity); V Q is the actual aircraft velocity; 
the error, £ , is the difference of the two; and is the 
Increment of elevator angle from an arbitrary reference. 

It was the desire of the authors of this investigation 
to obtain a mathematical relationship of the pilot's out- 
put U) versus the input (£ ) . This mathematical relation- 
ship,^, or transfer function, was desired in the 
simplest form that would give reasonably accurate results 
for a smooth air speed change. 

The authors felt that it would be possible to express 
the pilot's actions in a very simple linear relationship, 
&= a£ + b£ ; as the pilot sees the error and derivative 



of the error merely by watching his air speed Indicator 
reading and rate of change of reading. With practice, 
the pilot knows approximately where to position his 
elevator for a desired air speed change, and he positions 
it accordingly. 

Furthermore, there exists the pilot's Integrating 
ability, which physically corresponds to the pilot's 
desire to change air speed In a reasonable length of 
time. This ability Is also a factor which might effect 
the pilot's actions, and as such, should be considered. 

A number of basic assumptions were made in order to 
facilitate arriving at an expressable result. It was 
assumed that the pilot had only visual realization of the 
error or any simple function of the error, such as the 
rate of change, acceleration or integral of the error. 

No account was taken of stick force effect, as It was 
extremely small In the test aircraft; however, It Is felt 
thatthls effect definitely Influences pilot reation In 
extreme speed changes. The pilot's time delay was recognized 
to be of the order of one- third to three-fourths of a 
second. It is a variable quantity depending on the task 
that the pilot happens to be performing, and as such is 
extremely difficult to evaluate. No means were at the 
authors' disposal to record time delay; hence for the 
purpose of this investigation, It was ignored. Throughout 



this Investigation, the pilot was considered simply a 
component part of a servomechanism, or a "black box", 
and physiological relationships such as nerve response 
were not considered. 



PROCEDURE AND DISCUSSION 



The flight procedure for this phase of the Investiga- 
tion was Identical to that of the aircraft transfer function 
determination, with the Important exception that the air- 
craft speed was changed from approximately 130 mph true to 
110 mph true through the exclusive use of the elevator - 
power settings remaining constant. This speed change was 
made In a smooth, expeditious manner. Time histories of 
calibrated air speed and elevator angle were recorded for 
over 20 runs, two of which are Indicated In Figs. 10 and 
11 . 

The speed was changed on all runs In smooth air con- 
ditions at 5,000 feet mean pressure altitude. The pilot 
practiced before each series of recorded runs In order to 
make smooth speed changes void of overshoots or undershoots; 
hence the changes were not of an entirely random nature. 

The Initial effort to obtain a mathematical relation- 
ship between ^ and 6 was a numerical analysis of £, 6 and G 
in an effort to fit a curve of the order <3 = a 6 -h b<5 . 

The constants "a" and "b M were determined as indicated In 
the Sample Calculations. This equation necessitated 
assuming & to go Immediately to a finite value at t - 0. 

The Initial error was considered to be 20, and the final 
elevator angle and error was zero. 

These results were considered poor; so an electronic 
analog computer means, as Illustrated in Fig. 5 was devised 



for solving a second order differential equation In % and £ 

(aD 2 + bD + l)£ - (a'D 2 + b'D <-c')£ 

simultaneously with the aircraft equation as determined In 
Part I by varying the coefficients of the human equation 
to optimize each coefficient and arrive at a near perfect 
curve fit of i and £ for a reference air speed change of 
the test aircraft. 

With run 21?0 as a reference, these unknown coefficients 
were varied by adjusting the proper potentiometer settings 
and their Input resistances until the optimum curve fits In 
<f and £ for aircraft and computer curves were obtained. 

An original setting was made on each variable potentiometer 
to make both error and elevator curves stable. Following 
this, each variable potentiometer was adjusted until It 
gave optimum results; then the same procedure was repeated. 
The effect, c of varying each unknown constant was noted, and 
this appreciably assisted in obtaining a reasonably rapid 
curve agreement. The final result was checked with run 
2167. The constants "a" and "a' n were quite small and had 
little effect on improving the curves. 

In the aforementioned computer set-up, the original 
elevator angle and error were assumed to be non-zero and 
the final steady-state values zero. Initial conditions of 
+■ 20 mph error and +-. 65 ° elevator angle were set Into the 
computer problem. 



The optimized transfer function obtained was, 



i = H, » - (.01D 2 + + .0^2) 

6 (.01D 2 4 3 . 15 D + 1) 



and Ignoring the second derivative terms, as their 
effect Is negligible, yields, 




(■VSTO + .042) 
(3.15D + 1) 



The characteristic equation for the simultaneous equations, 
H f and A, was investigated and found to be, 

3o67d3 + 1.389D 2 t .53130 + .0775 - 0 



the roots of which are, 

A,--. .195; A 2>3 =- .0921 t .3i 7 i 

An attempt was made to analyze this result, but as 
the roots are of the same order of magnitude, a thorough 
analysis in terms of the constants of the transfer function 
was impossible. It is indicated, however, that the pilot 
in changing air speed is modifying the phugold mode by 
adding a comparatively highly damped convergent mode and 
increasing the phugold' s damping and frequency. 

The possibility of the pilot's integrating ability 
influencing the transfer function was next considered. 

This gives an equation of the form, 

ai+cf = a '<£ t- b'£ 4- c *J C dt 



The analog schematic for this determination Is shown 
in Fig. 12 with the final optimized potentiometer set- 
tings Indicated in Table II. The optimized equation 
obtained was. 



<T = 



6 



H, 



(.436D + .070?+/ 
(1.602D ■+■ 1) 




The characteristic equation for the simultaneous equations, 
H 2 and A, was investigated and found to be, . 

D U + .6? 8D3 +- .2981D 2 -h .0566D + .00279 = 0 
the roots of which are, 

\ = - .225; ^2 - - . 0 ? 4 ; ^ 3 _ 4 = - .188 t .3541 

♦ 

These roots indicate modes similar to the second 
order result, but contain an additional lightly damped 
convergent mode arising from the presence of the Integral 
term. 

In determining the aircraft transfer function for 
the integral solution, an initial velocity of +20 mph and 
a S of zero were assumed. This gave an equation of the 
form, 

<f + a £ + b £ c <£ - K 



For the error and its derivative to go to zero, the steady 



state perturbation elevator angle must equal K, or 
K. A value of S- - . 65 ° (<f for 20 mph change) 

was therefore set into the summation point for the 
aircraft transfer function. This in no way affected 
the dynamics or solution to the problem, as it merely 
shifted the arbitrary zero for & . 



RESULTS 



In analyzing the results of this phase of the invest- 
igation, note should be taken of the physical aspects of 
the results obtained. The optimized second order differ- 
ential obtained was: 

(.01D 2 4 - 3.I5D + 1 )/ = - ( 01 D 2 + .355D -h . 042)6 
This result indicates that the inertia effect of the 
elevator is small and can safely be neglected, and that 
even though the pilot is able to detect error acceleration, 
he changes his speed primarily according to his error and 
error derivative determination. The necessity of the 
inclusion of the & term is indicated due to physical 
limitations of the pilot to abruptly move his elevator 
control and his reluctance to move it too rapidly in 

Sfc 

consideration of aircraft stresses and pilot-passenger 

% 

comfort. The computer and aircraft reference error curves 
Indicated in Figs. 6 and 7 are excellent, but the elevator 
curves indicated in Figs. 8 and 9 do not agree too closely. 
They are within .05 degree and 1 second, however. 

By including the pilot's integrating ability the 
integro-dlf ferential equation obtained after optimizing 
the constants was \ 



lx. 00522) 

D 



(1.602B +1)S 



= -( .fc 36 D + .o?-,3f 



where "1/D" is the integral. The error curves of 
computer and aircraft reference were good, having a 
maximum variation of .3 mph, and the elevator curves 
were excellent, favorably comparing in both magnitude 
and phase. These curves are shown in Figs. 13 and 14. 

The algebraic sign of all terms of both the second 
order result and the integro-differential result should be 
considered. The error term in both human equations indi- 
cates to the pilot the proper direction to move his 
elevator (positive error, negative elevator angle). The 
error derivative acts as a limiter and tells the pilot 
he is approaching his desired air speed rapidly (for a 
high derivative) , and therefore to decrease his elevator 
angle (negative £, positive elevator angle. £ is always 
negative for a decrease in air speed). A positive 
integral of* the error tells the pilot not to linger in 
changing his speed, and also indicates the pilot's ability 
to detect how long the error has existed (positive integral, 
negative elevator). It is felt by the authors that the 
pilot not only has the ability to integrate but also the 
ability to reevaluate this integral from time to time. 



CONCLUSIONS AND RECOMMENDATIONS 



In arriving at an expressable pilot's transfer 
function, three attempts were made to obtain a result In 
the form of a differential or integro-differential equa- 
tion. The first, or simple system gave poor results and 
was not considered. The second order system results were 
good, and it must be concluded from these results that 
the second derivative terms are negligible. The Integro- 
differential result was excellent and reasonable, and it 
is concluded that with further study and investigation, 
the pilot's Integrating ability can further be explained. 

It is the authors' belief that the results of this 
investigation are quite limited in scope, as they were 
obtained for one pilot changing speed in one aircraft a 
definite increment under ideal conditions. In consider- 
lng a servomechanism system, however, and neglecting 
nonlinearities, test runs in different aircraft, with 
different aircraft transfer functions, should give close 
results to those obtained in this report for comparable 
air speed changes. 

A more complete Investigation should Include a 
frequency response or impulse integration analysis of 
various, random air speed changes, as outlined in Chapter 
11 of Ref. 8, to verify the pilot's response in terms of 
only an air speed error, its derivative and its integral. 



SAMPLE CALCULATIONS 



In order to fit a curve of the form, 

<£- &€ + b£ 

two simultaneous equations were formed from two sets 
of data points. Using known values of 6 % <S , and C 
one equation was formed from an averaging of the sum 
of even numbered seconds, and another from averaging 
odd numbered seconds. These two equations in Tt a" and 
n b rt were then solved simultaneously for a and b as indi- 
cated. 
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Dividing by five 
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+ 1.56 




42.3a 


- 9.476b 



or, +.316 = 3.460a - 1.09520b 



Determinant of coefficients, D, 



6.60 

0.46 



D - + 2.24761 

D a - - .1652104, a - - .0735005 
D b = -1.11220, b = - .494072 

Therefore, 

-.07350056 -.494072 



1.7442 

1.0952 



PART III 



CONSTRUCTION AND TESTING OF 
AN ELEVATOR CONTROL SIMULATOR 



THEORY 



The equations of aircraft longitudinal motion used 
in this section of the report are the simplified equations 
of Lift, Drag and Moment developed in Chapter 10 of Ref. 

1. These non-dimensional equations are: 
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where u.ol, and 9 are increments of velocity, angle of 
attack and pitch angle, and the operator "d" is d()/dt, 
where time, t, is in seconds. 

These three equations of motion are the stick-fixed 



equations, and as such do not take into consideration ele- 
vator hinge moment parameters and inertia effects of 
freeing the elevator. However, as the stick fixed phugold 
can be very closely approximated by the use of a positive 
control stop in the aircraft, and as the aircraft was 
small with light control surfaces it is believed that 
stick - free effects are thereby minimized in normal 
flight maneuvers. 

These equations of motion are in the three un- 
knowns, u,oC , and£* , the first of which is the non-dimen- 
sional output variable for the problems of this report. 



The 



quanti ty 



t 



x 2 - 



t 



is the forcing function or input 



variable for the problem. 

If these simplified equations of motion accurately 
represent the motion of an aircraft, then an elevator 
forcing function of a particular form should cause u,^, 
and 0 to vary In the same manner In both the actual 
aircraft and In the equations. This Is, of course, assum- 
ing that the aerodynamic and stability parameters used In 
the equations are correct. Also the converse should be 
true. If one variable, say u, were to be manipulated 
In a particular manner In the air, and the forcing 
function measured, then a similar manipulation of u In 
the equations of motion should reveal an Identical forcing 
function. 



PROCEDURE AND DISCUSSION 

The aircraft center of gravity and moment of Inertia 
were obtained by weighing the aircraft and oscillating It 
about its jack points. This procedure was carried out with 
the aircraft in its flight configuration - pilots positioned, 
fuel aboard, Instrumentation installed, and gear and flaps 
up. A schematic diagram, together with the calculations 
performed, may be found In the Sample Calculations section 
of the report. 

Stability and aerodynamic parameters required for 
the equations of motion were determined by calculation or 
from Refs. 3 or 4 as Indicated in Appendix B. For the 
purpose of these calculations, it was assumed that the air- 
craft was flying at 5000 feet pressure altitude at a true 
air speed of 110 mph. This was done in order to conform 
with the flight test conditions In Parts I and II of the 
report and to allow comparison with the data of those 
sections. These parameters when substituted into the 
simplified equations of aircraft motion shown in the theory 
section of Part III give: 

LIFT : .400 u+(.213 + d)^ - d-©- = 0 

DRAG : ( .038 +d)u - .0420-^ ^.200^ = 0 



MOMENT : ( ,2?5 +-. 450 d)*c + ,io d-0- - .0642d 2 -e- = .^55 



These three equations were set up on an electronic analog 
computer as indicated In Pig. 18 and solved simultaneously 
using a step function input voltage corresponding to an 
Incremental step function of .8 degrees elevator angle. 

The velocity response to this step function input 
was recorded on a Brush Pen Recorder. It was desired to 
record the actual velocity change, V, Instead of the 
non-dimensional incremental velocity, u, so the relation- 
ship AV • Vu was used to calibrate the velocity trace. 

As a step elevator input of .8° was used as a forcing 
function the velocity response of the computer system 
could be compared with the actual response of the aircraft 
during Run 2172 as shown in Fig. 3 and discussed in Part I. 
The actual aircraft had a period of 36 seconds and a time 
to damp to half amplitude of 25 seconds during this run 
while the solution of the simplified equations of motion 
by the computer indicated a period of 31 seconds and a 
time to one half amplitude of 58 seconds. These discrep- 
ancies may be attributed to the following reasons: the 

aircraft is actually non-linear in the speed range con- 
sidered, that stick fixed conditions did not actually 
exist due to control cable stretch or flexibility, and 
that propeller effects were not included in the computa- 
tions for the Cp used in the simplified equations. 

Ref. 5 indicates the necessity of including this last 
factor. 



In an effort to make the computer velocity response 
to a step elevator Input more nearly correspond with 

d 

that of the aircraft the value of-~- was changed to .1 
(an unreasonable value). This gave a period of 31 seconds 
and a time to damp of 29 seconds which was a close enough 
approximation for the purposes of this section. The 
expression for the aircraft transfer function as obtained 
in Part I of the report was not used, as it was desired 
to check the classical equations as indicated above. 

A mock-up to simulate the elevator control of the 
test aircraft was constructed and is illustrated in Figs. 
15 and 16. The mock-up consisted of a wooden frame 
supporting a Navion control column. The fore and aft 
movement of the control column was restrained by 4 elastic 
strips which roughly simulated the stick force required 
to move the elevator in flight. The control column was 
connected through a pulley system to a precision poten- 
tiometer which enabled picking off positive or negative 
voltages corresponding to fore and aft control (or down 
and up elevator) movement, respectively. The ratio of 
the control displacement to elevator angle change in the 
mock-up was adjusted (as explained later) to approximately 
the same value as in the actual aircraft. A meter whose 
deflection was proportional to u (and thus to^V) was 
displayed in front of the pilot immediately above the 
control column on the mock-up. One side of this meter 



was connected to ground on the computer and the other 
side was connected to the u output of the computer. 

This meter is shown in Figs. 15 and 16 and in the schematic 
Fig. 17. It should be noted that the scale of the meter 
read from -100 to-t-100 with 0 in the center position. Thus, 
deflections moving toward the right indicated increases in 
air speed while those moving to the left indicate decreases 
in air speed. When the stick was in its neutral position 
there was no S input and consequently no u output from 
the compute^ and the u meter read zero. 

In order to record the time history of the change in 
velocity and the pilot's movement of the control column, 
a two channel Brush Pen Recorder was connected to the 
computer to record the input b and the output u*LV. 

Sc 

The voltage supply across the potentiometer was adjusted 

% 

so that a control movement which produced 1° of elevator 
angle change In the airplane caused approximately 1 volt * 

1 degree input to the computer. The controls of the 
recorder were adjusted to give the scale factors indicated 
in Figs. 19 and 20. Thus a 60 unit deflection on the u 
meter corresponded approximately to a 13 mph change in 
velocity while a 120 unit deflection of the meter was 
roughly a 26 mph change. It should be noted that these 
calibrations are close approximations and not exact. 

The pilots practiced "flying" the mock-up by starting 
with the control column in the neutral position (u meter 



reading zero) and attempting to decrease the meter read- 
ing (and thus the velocity), through movement of the 
controls, to the -60 position ( k V = 13 mph) in a smooth 
and expeditious manner and arriving with no overshoot 
or undershoot. A recorder tra.ce of a run of this type 
is shown as Pig. 19. Other runs were performed in which 
the pilot started from a stabilized value of meter de- 
flection of + 60 and attempted to decrease it to -60. 

This corresponded to a decrease In air speed of 26 mph. 

A record of this type air speed change is shown in Pig. 

20. In all runs it was found very difficult to move 
the controls smoothly and evenly due to the large amount 
of stlctlon in the mock-up. It was never possible to 
arrive at the final meter reading without some oscilla- 
tions . 

Observation of one pilot "flying" the computer by 
another, plus discussion between pilots as to how they 
flew the computer yields the following information. The 
initial control movement made in both types of runs to 
Initiate the speed change was one which, in the memory or 
experience of the pilot, corresponded to that which he 
had made in the test aircraft to produce a similar air 
speed change. Control movements subsequent to the initial 
movements varied in each type of run. In the case of 
the 60 unit change of the u meter the initial stick 



position was held until about 50 units of change had 
occurred, at this time about half the control displace- 
ment was removed in order to avoid overshoot. This was 
almost immediately followed by a backward movement of the 
control column to avoid undershoot. Then followed a 
random motion of the control in an attempt to quiet the 
oscillations of u about the final value. In the case of 
the 120 unit change the initial control movement was held 
until the needle slowed perceptibly (which occured short 
of the desired change), at this time an additional rear- 
ward movement of the controls was made to further increase 
the value of u. This was followed by random motion of 
the elevator to stop the oscillations about the final 
value. 

* 

An examination of the records of these two type 
runs. Figs. 19 and 20, immediately reveals the following 
things. The pilot handled the elevator control in the 
mock-up far more roughly than he did in the actual air- 
plane. This can probably be explained by the lack of 
"G" resulting from the control movement, the failure to 
properly simulate the control forces, and the large amount 
of static friction in the system. The pilot reached his 
final meter reading within 13 to 15 seconds after his 
initial control movement regardles of whether he was 
changing air speed 13 or 26 miles per hour. The pitch 
rate of the aircraft would be higher in making large 



speed changes than in the smaller change. The pilot 
experienced difficulty in settling down to his new meter 
reading without oscillations. 

Before attempting to explain these phenomena it is 
necessary to consider the velocity response to an elevator 
step function in both the test aircraft and the computer 
equations in more detail. As the equations of motion in 
the computer are linear, an elevator step function input 
will cause the phugold mode to be excited and the 
amplitude of the oscillations will be in a definite ratio 
to the size of the step input. If a bigger step is applied, 
a bigger oscillation will result, but this ratio will 
remain constant, and, in addition, the period and the time 
to damp to half amplitude will remain the same. This 
is illustrated below. 
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The aircraft is not linear over the range of speeds tested 
in Part I, however, the same relationship is approximately 
true as is shown in Table I and Figs. 2 and 3* 

As mentioned before, the aircraft equations of 
motion used in the computer did not exactly match the 
actual aircraft transfer function found in Part I. Con- 
sequently, Figs. 2 and 3 can be used only for qualitative 
comparison with the mock-up-computer system results. 
However, it can be seen that in "flying" the computer, 
the pilot put on an initial step function and waited to 
see how it affected the velocity response. The effect he 
witnessed was essentially the first oscillation of the 
phugoid mode. In the case of the 13 mph decrease in air 

speed (Fig. 19), he generally put on an Initial step 

. * 

elevator deflection that caused a phugoid oscillation of 
too great an amplitude which would have carried him 
beyond his desired speed had he not applied opposite 
control as he realized he was approaching his final speed 
too rapidly. The opposite control, which was applied at 
about 12 or 13 seconds, together with the decrease in 
the rate of speed change which occurred as the phugoid 
oscillation approached its maximum (at about 15 or 1 6 
seconds) enabled him to reach an airspeed close to his 
desired one. However, at this time the phugoid oscilla- 
tion caused the air speed to commence increasing. This 



necessitated a large amount of up elevator (back 
control) to stop. Thereafter, the pilot's control 
movements were "desperation" movements to stop the 
succeeding small oscillations. In the case of the large 
air speed change (Fig. 20), the pilot usually applied an 
initial step Input that was not sufficient to cause a 
phugold oscillation which would take him to his desired 
air speed. Thus, as he watched his air speed needle 
move, the rate of change was satisfactory for the first 
12 or 13 seconds until the rate of change of air speed 
began to decrease due to the phugold approaching the 
maximum of its first oscillation, and it seemed that he 
would not reach his desired speed. As a result, he 

added more up elevator which brought him to his final 

* 

speed and also partially compensated for the change in 
phugold oscillation tending to increase his air speed. 

As this increase in speed due to the phugold progressed 
he added more up elevator. By this time he had the mock- 
up fairly well under control, and his subsequent movements 
were not as radical as in the smaller air speed change. 

From the preceedlng discussion it may be seen that 
the primary reason the pilot reached his new air speed 
in 13 to 15 seconds was due to his step function elevator 
input. This input Initiated a phugold mode whose first 
oscillation took 15 to 16 seconds to reach a maximum, 
and whose amplitude was approximately that of the desired 
air speed change. The oscillations of the air speed about 



the final desired value, end the associated radical 
control movements occurred largely as a result of the 
reversal of the phugold oscillation which tended to 
cause a reversal In the direction of air speed change. 

The Increased pitch rate In the case of the larger 
air speed change may be considered to stem from two causes. 
First, the pilot when flying the computer had no sensation 
of "G" force and consequently had no hesitation in apply- 
ing an elevator step change of considerable magnitude 
immediately, rather than applying It gradually as he did 
1: ' he aircraft. Secondly, the pilot was reluctant to 

put on as much control deflection in the case of the 
small speed change as he would use In a larger change for 
fear he would overshoot. 

*3r 

An attempt was made to correlate the time history 
*• 

of velocity and elevator response as recorded from the 
mock-up runs, Figs. 19 and 20, with those recorded during 
the smooth air speed changes In the test aircraft in 
flight (Figs. 5 through 11). Because of the difference 
In the values of air speed changes and elevator Inputs 
on the mock-up and In actual flight this was difficult to 
do. However, It was noted that the velocity response curves 
for runs 216? and 2170 were almost identical In character 
with the 26 mph mock-up changes during the first 9 seconds. 
After this time, the flight test velocity response 
commenced to level out as the pilot approached a final air 



speed 5 ®ph less than that made on the mock-up. This 
Initial coincidence does show the similarity of the 
Initial transient behavior of the phugold in each case. 

No comparison between the elevator movement could be 
made as the movement of the elevator controls in the 
mock-up was far more erratic and rough than when in the 
air. Had the mock-up been more carefully constructed so 
that the control forces were more nearly correct, and had 
the static friction been removed from the system, smoother 
elevator control would have been realized. 

When flying the aircraft and making smooth air 
speed changes, the pilot did his best to concentrate on 
the air speed Indicator in an attempt to have this as his 

only source of information when making the air speed change. 

* 

He did, in fact, have other means at his disposal to 
help evaluate the manner in which he was performing the 
maneuver. The 11 seat of his pants" gave him information 
as to whether or not the "G"load on the aircraft was 
excessive or uncomfortable and his peripheral vision 
Informed him of the aircraft attitude and pitch rate. 

With these added sources of information at his disposal 
he was able to overcome changes in sign of the phugoid 
oscillation with no conscious effort on his part by merely 
exerting pressures on the controls rather than making 
actual control displacements. It is felt that if the pilot 
had had additional information concerning pitch rate 



at his disposal while making the mock-up air speed changes 
his control positioning would have been less erratic,, 

Another, and still better, method of making a more 
reasonable correlation between the flight test results 
and the computer mock-up runs is the following. Conduct 
the flight tests with the pilot essentially flying blind, 
and with all instruments covered except the air speed 
indicator. This would remove all information sources except 
"G" and that portrayed by the air speed indicator. Instead 
of the classical equations of motion being used in the 
computer, they should be replaced with the actual aircraft 
transfer function as determined in Part I. If these 
techniques were followed and a mock-up of better design 
were incorporated, it is felt that a reasonable and inter- 
esting comparison might be made. 



BESULTS 



The analog computer solution of the classical 
equations of stick fixed motion were found not to accurately 
represent the aircraft velocity response when the speed 
change is of the order of 10 to 25 mph. This may be 
attributed to the following: the aircraft is not linear 

in the speed range considered; that stick fixed conditions 
did not actually exist due to control cable stretch or 
control flexibility; that propeller effects were not con- 
sidered in the computation of the drag coefficient. The 
phugoid period as computed by the classical ecuations was 
within 15 % of the value as determined by flight test. The 
time to damp to half amplitude computed by the classical 
method was greater than that found in flirht test by a factor 
of 2. These values of period and time to damp could be 

a 

changed to more nearly agree with the actual aircraft by 
changing the value of drag coefficient by a factor of 2 
to 3. In future work of this sort it is felt that the 
value of the aircraft transfer function as determined by 
flight test in Pa.rt I would be more satisfactory than trying 
to adapt the classical equations of motion. 

It was found that insufficient attention had been 
devoted to the mock-up construction. The static friction 
in the control column was too high. This could be reduced 
by the use of proper bearing surfaces. The control forces 
were not similar to those encountered in the actual air- 
craft. Better spring restraints on the control column 



and a type of "Q box" would remedy this situation. The 
above considerations added materially to the difficulty 
of "flying the computer. The method of presentation of 
the air speed data on the mock-up was completely adequate, 
and similar to that which the pilot encountered in normal 
flight. 

Changing air speed on the mock-up smoothly and with no 
overshoot or undershoot proved extremely difficult. In 
addition to the mock-up limitations enumerated before, 
certain other factors contributed to this difficutly. The 
pilot positioned his controls initially as a step function 
elevator input in accordance with his memory as to how he 
had displaced the controls in the aircraft. Then, in 
the absence of any information as to "G" or pitch rate, 
he was forced to adopt a "wait and see" attitude as the 
phugoid raod-e. progressed. In general, the amplitude of the 
first oscillation of the phugoid mode which was excited by 
the initial control deflection was such as to reasonably 
approximate the final air speed change desired. Subsequent 
oscillations about the final value of air speed seemed to 
be closely associated with the phugoid oscillation revers- 
ing and causing an increase in air speed. Control movement 
in this area seemed to be random in nature. Pilot opinion 
indicated that it was easier to stabilize the mock-up 
after a large air speed change. This may be due to the 
fact that the pilot had more experience in making large 
air speed changes on the computer, or that, as his initial 



control Input In the case of the large air speed change 
was not sufficient to cause a phugold oscillation that 
would overshoot his final airspeed, he was forced to add 
additional up elevator to get to his final speed and this 
additional elevator helped overcome the tendency for the 
air speed to Increase as the phugold oscillation reversed. 

As the motion of the elevator in the computer mock- 
up runs was so erratic, no attempt was made to compare 
them with the elevator motion in flight test. It was felt 
that with an Improved mock-up design, and a flight test 
procedure that limited the sources of pilot information to 
the air speed indicator and to "G n would allow intelligent 
correlation of the data. 



CONCLUSIONS 



The classical equations of stick fixed, aircraft 
motion do not accurately represent the aircraft velocity 
response when the speed range Is of the order of 10 to 
25 mph. 

It would have been better to have used the aircraft 
transfer function as obtained in Part I in lieu of attemp- 
ting to adapt the classical equations. 

Static friction in the mock-up control was too high. 

Control forces in the mock-up should have more nearly 
simulated those of the test aircraft. 

In attempting to change air speed on the mock-up the 
pilot applied a step function elevator input which correspond 

ed in his mind to the stick displacement required for a 

* 

similar air speed change in the test aircraft. This elevator 
deflection excited the phugold mode, and the magnitude of 
the first oscillation brought him very close to his desired 
air speed change. 

Oscillation about his final air speed with the 
accompanying erratic control movements represent an effort 
to nullify the effect of the continuing phugold mode. 

No comparison of the elevator movement could be made 
between the in-flight air speed change and the mock-up 
speed changes. 

No comparison of the elevator movement could be made 
between the in-flight air speed change and the mock-up 
speed changes. 



The absence of other information such as acceleration 



and pitch rate made the pilot's control movements rough 
and erratic, and caused him to adopt initially a "wait and 
see" attitude followed by a slam bang "I can stop that" 
philosophy. 

In any future investigations it is recommended that 
the mock-up construction be Improved, and that either more 
information be displayed on the mock-up or that flight 
tests be conducted where .the pilot is flying blind with 
all i:;struments except' the air speed indicator covered in 
order to allow for more intelligent correlation of air 
speed and elevator data. The aircraft transfer function 
as determined in Part I should be used in the mock-up 
computer system in future tests. 




DETERMINATION OF AIRCRAFT CENTER OF GRAVITY 



AND MOMENT OF INERTIA 



Veight on jack points: 2,580 lbs. 

weight on tail skid. : 168 lbs^ 

Total \ ei£it : 2 , 746 lbs . 

(2,5^0) (89) + . (168) (315*32) - 2,748X 

It . 229,900.+ g2,?0Q . 102.9 inches 

2,748 

$ K»C - 102, 6 3 ~35 ? * 66 - 28.3S5 MAC , X<. g - .281 



Assuming an undamped oscillation, 



P^kl 2 

J ack " UK 2 



where "k" is the 



spring constant, n l u the spring lever arm and rt P n 
the period of the oscillation. 

Test data: k « 70.65 lb./in.; P av = 



.666 sec 



